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Abstract 

We address Standard Model predictions for flavor-changing radiative 
transitions of the pseudoscalar charm mesons. Short-distance con- 
tributions in D radiative transitions are contrasted with those in B 
decays. A full analysis is presented of the c — ► u + 7 electromagnetic 
penguin amplitude with QCD radiative corrections included. Given 
the importance of long-range effects for the charm sector, special at- 
tention is paid to such contributions as the vector dominance and pole 
amplitudes. A number of two-body final states in exclusive charm ra- 
diative decays is considered and the corresponding branching ratio 
predictions are given. 



1 Introduction 



Important milestones in the study of the 5-quark system were reached with 
the recent observations of both the exclusive decay B — > K* 

B B ^ K * 7 = (4.5 ± 1.5 ±0.9) x 1(T 5 , (1) 

and of the inclusive transitionB 

B 6 ^ S7 = (2.32 ±0.57 ±0.35) x 10~ 4 . (2) 

To first approximation, these flavor- changing radiative decays can be inter- 
preted at the quark level in terms of the b —+ s'j transition. The Standard 
Model allows for such a process by means of a one-loop penguin-type ampli- 
tude. Within errors, agreement of the measured branching ratios and Stan- 
dard Model predictions appears to be reasonable. The small magnitudes of 
these branching ratios indicate just how sensitive experimental probes of b- 
quark hadrons have become. However, an outstanding question is the size of 
the long distance contributions 0, || to such radiative B decays relative to 
the short distance penguin amplitude. This issue must be addressed in order 
to establish the viability of determining the value of the ratio of CKM matrix 
elements, | Vtd\/\ Vt s \, from a measurement of the ratio of exclusive branching 
fractions Bb^p^/Bb^k*^- 

Important as they are, the above measurements by no means exhaust 
the set of interesting problems. It has become increasingly evident that 
the database for charm hadrons is also in a state of rapid expansion, and 
that physically important levels of sensitivity are being achieved. Perhaps 
the most impressive example of this to-date is the recent observation of the 
nonleptonic decay D° — > K + n~, with branching ratio0 

,L> "- Kir - 0.0077 ±0.0025 ±0.0025 . (3) 



B D °- 



This transition has been interpreted as evidence of a doubly Cabibbo sup- 
pressed transition rather than of D°-D° mixing. 

The discussion in this paper will be directed towards a somewhat different 
aspect of charm physics, the flavor-changing radiative decays. These transi- 
tions require thejoint occurrence of weak and electromagnetic interactions. 



3 tne ion 
1§,0,§ 



From Table lQ'liJ'ltJ, we see that no such events (involving emission of real 
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or virtual photons) have yet been observed. However, these decays are an 
active area of study and data gathered in ongoing fixed-target experiments 
are establishing markedly improved bounds. Our objective in the analysis 
to follow will be to provide up-to-date predictions for flavor-changing ra- 
diative transitions of charm systems. Since the experimental situation for 
charm mesons is at present more favorable than for charm baryons, we shall 
restrict our attention to the former. Even with this restriction, it is a tall 
order to supply accurate theoretical values. It has become evident over a 
long period of time that theoretical calculations of D-meson weak decays are 
not particularly trustworthy, due in part to the absence of a rapidly conver- 
gent approximation scheme and also to the presence of significant hadron 
dynamical effects in the D meson mass region. Despite this, we feel that one 
can make some definite statements, such as the relative importance of long- 
range and short-range effects and of the various types of final states which can 
reasonably be anticipated. We shall base our analysis on a variety of theoret- 
ical techniques, from operator-product expansion and renormalization-group 
methods to more phenomenological approaches like vector-meson-dominance 
(VMD). Measurement of radiative charm decays would probe the long dis- 
tance contributions and thus provide further insight in the extrapolation of 
calculational techniques to the B sector. 



Table 1 Status of Electroweak-induced Charm Decays 



Mode 




Branching Ratio 


D° — »■ 


A 


< 1.4 x 1(T 4 


D° — > 




< 2.0 x 1CT 4 


D° — > 


K°e + e- 


< 1.7 x 1CT 3 


D° — > 


p°e + e~ 


< 4.5 x 10~ 4 


D° — »■ 




< 8.1 x 1(T 4 


D+ -» 


• 7r + e + e~ 


< 2.5 x 1(T 3 


D+ -» 




< 2.9 x 1(T 3 


D+ -» 


■ K + e+e- 


< 4.8 x 1(T 3 


D + -» 


■ K+fi+fj,- 


< 9.2 x 1(T 3 



Let us summarize the contents to follow. In Section 2, we consider the 
short-range component in radiative charm decays, primarily the charm coun- 
terpart of the penguin amplitude which dominates the radiative 5-meson 
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decays. In addition to addressing c-quark physics, our analysis contains a 
purely theoretical advance by removing an unnecessary assumption made in 
earlier studies involving 6-quark applications. Section 3 begins our analysis 
of the so-called 'long-distance' contributions with an analysis of pole dia- 
grams, which are induced by the weak mixing of pseudoscalar and/or vector 
charm mesons with noncharm states. In Section 4, we continue our study 
of long distance effects by turning our attention to a study of VMD ampli- 
tudes. Our conclusions and recommendations for future studies are given in 
Section 5. There is also an Appendix in which the applicability of VMD to 
certain light-meson decays is commented on. 

2 Short Distance Contributions 

Examples of diagrams which mediate the short- distance transition ampli- 
tudes for radiative charm decay are depicted in Fig. 1. They have in com- 
mon that the photon emission occurs in a region of spacetime determined 
by the propagator of the W-boson. In view of the large W-m&ss -Mw, this 
region has a very limited extent compared to the length scale of the strong 
interactions, hence the name 'short- distance'. Looking ahead, our conclusion 
regarding such short-distance amplitudes will be that in radiative decays of 
charm mesons they are small relative to long-distance effects, even though 
they receive large enhancements from QCD corrections. As described earlier, 
this is of course in stark contrast to B decay. 



Figure 1. Short- distance Effects 

In recognition of the importance attached to the electromagnetic penguin 
transition in radiative B decays, we give a brief pedagogical comparison 
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between the role played by this effect for B and for D decay. The two 
transitions in question are given at the quark-level by 

6(p,A) - S (p',A')+7(q,^) 

c(p,A) - u(p',A') + 7(q^) • (4) 

To highlight the crucial role played by the quark masses and CKM matrix 
elements, let us at first ignore the effect of QCD radiative corrections. The 
relevant Feynman diagrams are then depicted in Figs. 1(a), 1(c) and the pen- 
guin amplitude for the transition of a heavy quark Q to a much lighter quark 
q and an on-shell photon is given by® 



eG F 



\F 2 {xi) u q (p', A')e Mt (g, \)a^q u [m Q P R + m q P L ]u Q (p, A) 



4V27T 2 

where Pr(l) are the right (left) -handed helicity projection operators, Xi = 
mj/Myf, Aj = V; s *Vib for Ab^ sl and \ = V c *V ui for A C -^ UJ . The function F 2 
gives the contribution of each internal quark to the electromagnetic penguin 
loop, 



F 2 (x) = 

x 3 — 5x 2 — 2x 3x 2 In x 



Q 



+ 



2x 3 + 5x 2 — x 3x 3 In x 
+ 4(x- l) 3 2(x- l) 4 ' ( ' 



4(x - l) 3 2(x - 1 

with Q being the charge of the internal quark. For b —>■ sj the sum is carried 
out over the quarks u, c, t and the term proportional to the s-quark mass in 
Eq. (H) is generally neglected, whereas for c — > uy, one sums over the quarks 
d, s, b and ignores the corresponding term proportional to the w-quark mass. 

Let us get acquainted with some of the numerical values. In Table 2, 
we first display the magnitude of the function F 2 and then fold in the CKM 
dependence for the b sy transition (we take m u = 5 MeV, m c = 1.5 GeV, 
mt = 174 GeV, and the central values of the CKM matrix elements as given 
in Ref. §). 

Table 2 Contributions to b —>■ s + 7 



Quark 


F 2 




u 


2.27 x 1(T 9 


1.29 x 10" 12 


c 


2.03 x 1(T 4 


7.34 x 10~ 6 


t 


0.39 


1.56 x 10~2 
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Dominance of the t-quark intermediate state is evident, even upon including 
the CKM factors. Its effect is so large that the other intermediate states are 
numerically negligible and hence are typically omitted. The corresponding 
situation is given for c — > wy in Table 3 (with m d = 11 MeV, m s = 150 MeV, 
and m b = 4.9 GeV). 



Table 3 Contributions to c — > u + 7 



Quark 


F 2 


\V c *V ni \F 2 


d 


1.57 x 1(T 9 


3.36 x 1(T 10 


s 


2.92 x 1(T 7 


6.26 x 1(T 8 


b 


3.31 x 1(T 4 


3.17 x 10" 8 



The amplitude for c — > w~) differs from that of b — > 57 in two important 
respects, (i) there is no single intermediate state which dominates, and (ii) 
the overall magnitude is much smaller. 

Neglecting the final state fermion mass, the QCD uncorrected decay rate 



^oL qi is given by 



aG% * 
128tt 4 Q 



(7) 



To obtain the branching fraction, the inclusive rate is scaled to that of the 
semi-leptonic decay Q — > q'£u. This procedure removes uncertainties in the 
calculation due to the overall factor of itlq which appears in both expressions, 
and reduces the ambiguities involved with the imprecisely determined CKM 
factors. Taking the above numerical values for the internal quark masses, and 
using the values of the semi-leptonic branching ratios as given in Ref . || , this 
yields 



+S7 



B, 



3a 


|^F 2 (x t )| 2 




g(m c /m b ) + - 


v ub 
v cb 


2 

Ygimu/rrih) 


= 1.29 xl(T 4 , 






3a \V; s V us F 2 (x s ) + V; b V ub F 2 (x b )\ 2 


271 ivy* 


g(m s /m c ) + - 


v cd 


2 

^g(m d /m c ) 



B 



B-,Xlv 



(8) 



B 



D+^X£+u , 
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= 1.39 x 10 

Here, the function g(x) is the usual phase space factor in semi-leptonic me- 
son decay, where constituent values of the final-state quark masses have been 
used. 10 The QCD uncorrected c — > wy transition is seen to have an unob- 
servably small branching fraction. 

We next examine the impact of the QCD radiative corrections on the 
above branching ratios. We begin by reviewing the calculation for the b — > 57 
transition, which will serve as the foundation of our subsequent discussion of 



c — ► uj. The QCD corrections are calculatedGJ via an operator product 
expansion based on the effective hamiltonian 

^-^E^W), (9) 

where the {Ok} are a complete set of renormalized dimension-six operators 
involving light fields which govern the b — > s transitions. They consist of 
two current-current operators Oi,2, four strong penguin operators O3-Q, and 
electro- and chromo- magnetic dipole operators O7 and 8 , 



1 = 

2 = 

3 = 


(c a -f^P L bp)(Sf3YP L C a ) , 
{CalfiPMisp^PLCp) , 

q 


4 = 


1 


5 = 


(s a -i^P L b a ) ^2(q/rfPRqfi) , 
q 


o 6 = 


q 


7 = 
Os = 


—^m b {s a a IJiV P R b a )F IJ ' v , 
lo7r 

-^-m b {s a a^T^P R b p )G a ^ 



(10) 



The above effective hamiltonian is then evolved from the electroweak scale 
down to n ~ nib by the Renormalization Group Equations (RGE). 

In the RG analysis, the Wilson coefficients are to be evaluated perturba- 
tively at the W scale where the matching conditions are imposed and then 
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evolved down to the renormalization scale /i. The expressions for the {ck} 
at the W scale are 



Ci, 3 _ 6 (M W ) 

c 7 (M w ) 



, 

-\F 2 {x t ) , 



c 2 (M w ] 
cs(M w ) 



-\D{x t ) 



with 



D(x) 



x 3 — 5x 2 — 2x 3x 2 lnx 
+ 



4(x- l) 3 2(x- l) 4 ' 
The solution to the RGE at the leading logarithmic order is given by 

c e k ff (fi) = U 5 ke (fi,M w )c £ (M w ), 



(12) 



(13) 



where ET^ denotes the evolution matrix in a five-flavor context and is deter- 



mined by 



U 5 (m 1 ,m 2 )kn = O ke 



(14) 

In the above we define rj = a s (m 2 ) / a s (mi) and di = 7^/2/? (not summed on 
£), where (3q — 11 — 2n//3 and 7 D = O -1 ■y( e ff^ T O is the diagonalized form 
of the 8x8 anomalous dimension matrix. We use the scheme-independent 
form of the matrix 7 6 -^, which is given explicitly in Ref. [HJ in terms of the 
number of Q = +2/3 and Q = —1/3 quarks present in the effective theory. 

Scaling again to the semi-leptonic decay, the branching fraction is now 
given by 



6a 

71 



v tb v; s 2 



v. 



cb 



I eff 
\ c 7 



g(m c /m b ) + -jj^g(m u /m b ) 



•B 



B^Xtv 



(15) 



cb\ 



The numerical values of the separate contributions to (//) are, with m t 
174 GeV (for illustration purposes), /i = m b = 4.87 GeV, and a s (M z ) 
0.124 as determined by LEP0, 



(//) = 0.670 c 7 (Mvk) + 0.091 ^(Mvk) -0.172 c 2 (M w ) 
= 0.670(-0. 195) + 0.091(-0.097) - 0.172 . 



(16) 



Taking the overall CKM factor in the branching fraction to be unity, and 
\V ub \/\V cb \ 



0.08, this procedure yields 
B^ S7 = (2.921°; 



x 10" 



(17) 



7 



The central value corresponds to fi — nib, while the upper and lower er- 
rors represent the deviation due to assuming /i = nib/ 2 and u, = 2nib, re- 
spectively. We see that this value compares favorably to the recent CLEO 
measurementB of the inclusive rate cited earlier in Eq. (|2|). When com- 
pared with the uncorrected result of Eq. @, the QCD corrections are seen 
to increase the branching ratio by roughly a factor of 2. 

We take this opportunity to reflect further on the size of the QCD correc- 
tions. Earlier estimates^ of these corrections found that the enhancements 
to the b — > S7 branching fraction were more than an order of magnitude for 
m t < Mw- This is because the effect of the QCD radiative correction to the 
weak vertex is to replace the GIM power suppression in Eq. (|5|) by a logarith- 
mic suppression. We explicitly illustrate this effect in Fig. 2, where we show 
the dependence of the C7 Wilson coefficient on the mass of a single internal 
quark using the calculational procedure described above. In the lower of the 
two curves, the dependence of c 7 determined at scale /1 = m w is displayed, 
while the upper curve corresponds to the evolved C7 evaluated at \i = nib- We 
see that cy(/z = nib) is a reasonably flat function of the intermediate quark 
mass, and that the corrections are substantial for light internal quarks, with 
an increase of 3-4 orders of magnitude in the rate for m q = 5 — 10 GeV. For 
the case of one heavy internal quark, e.g., b — > 57 with m t > Mw, we see 
that the GIM mechanism no longer plays such a crucial role and the QCD 
enhancements are not as dramatic. 
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Figure 2. Dependence of c 7 on Intermediate-quark Mass 

We now consider the case of radiative charm transitions. The |Ac| = 1 
effective hamiltonian can be written as 

HlP =1 = ^\ b j:cMOM, (18) 

with Xi = V^Vui as defined previously. The CKM structure of the operators 
differs dramatically from the b — > case. Here, both 0\ and O2 have two 
contributions which have approximately equal CKM weighting since X s ~ 
Xd- We stress that extreme caution must be exercised in order to correctly 
incorporate these terms. To be precise we explicitly separate 0\ and 2 into 
two operators according to their CKM structure, 

Ola = (u a ^ij,PLSp)(s^P L C a ) , lb = (u a ^P L dp)(d^P L C a ) , 

02a = (u a 'J^P L S a )(s l3 'J ,1 P L Cp) , 2 b = {UocltiPLda^dpYPLCp) , 

(19) 

and write the remaining |Ac| = 1 operators in a form analogous to their 
|A6| = 1 counterparts, 

3 = {u a 7 lx P L c a )Y^{qi3l' x PLqp) , 
q 
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4 = 


{Uofl^PLCp) 


<? 


o 5 = 




<? 


Oe = 






7 = 


1R 2 m c (u a 

1D7T 2 




Os = 


9s 

1C 9 m c (u a 
167r^ 


a, v T a ap P R c p )G^ 



(20) 



where the terms proportional to m u in O7 8 have again been neglected. Since 
the quantity —At has been factorized in Eq. ( |18|) above, the values of the 
corresponding Wilson coefficients at the matching scale are now 



Ci a (M W ) 
C2a(M W ) 







-A s / Aft 



ci 6 (M w ) 
c 2b {M w ) 







■Ad/ At 



(21) 



The values of the Wilson coefficients for cs-q(Mw) are the same as in Eq. (|TT|), 
and the coefficients cj^(Mw) are modified to 



c 7 (M w ) 
c 8 {M w ) 



^F 2 (x s ) + F 2 (x b ) 
^D(x s ) + D(x b ) 



(22) 



with each containing intermediate s-quark and 6-quark contributions. Due to 
the CKM dependence, ct$(Mw) now contain both real and imaginary terms 
which in principle must be evolved separately. We note that the real parts of 
the avdependent terms are numerically the same order of magnitude as the 
x b terms. Now we evolve the effective theory down to the scale [i ~ m c . This 
takes place in two successive steps; first, we go from the electroweak scale 
down to m b working in an effective 5 flavor theory, and then to /i < m b in 
an effective 4 flavor theory. Thi s pr ocedure is similar to what is performed 
for the I As I = 1 kaon transitions^ , where the effective theory is evolved to 
fi ~ 1 GeV in 3 successive steps. We then have 



2hic k ff {fi) 



U&n, m b )Ul(m b , M w )-Rec n (M w ) 
U^i/J,, m b )Uf n (m b , M w )Jmc n (M w ) 



(23) 
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and 

\c e 7 ff (fx)\ 2 = \-Rec e 7 ff {fi)\ 2 + \Zn4 ff {v)\ 2 . (24) 

The renormalization group evolution matrices U A and U 5 are determined 
as in Eq. ( fRD now using a 10 x 10 anomalous dimension matrix 7 e -^. We 
take the anomalous dimensions of the split operators 0\ a ~b and 0<i a -b to be 
exactly those for 0\ and O2, respectively, as the anomalous dimensions do 
not depend on the CKM structure of the operator. We use the form of 



as given in Ref . [13[] , taking care to keep n/ = 4 and 5 as needed. The relative 
numerical values of the contributions to C7 (/^) with fi = m c = 1.5 GeV are 

C7 (m c ) — 72ec7 (m c ) , 
= 0.458 c 7 (M w ) + 0.125 c 8 (M w )-0.312[c 2a (M w )+c 26 (M w )], 

0. I58f-G.2 I I x Kr 6 ) + 0.125f— 0.139 x 10~ 5a n Q1 ° ' ~ As ~ Ad 



= 0.458(-0.241 x 10~ 6 ) + 0.125(-0. 139 x 10~ 5 ) -0.312, (25) 

where the CKM unitarity condition A s + A^ = — has been used to simplify 
the final term. Incidentally, it should be stressed that the choice of —\b as a 
prefactor in Eq. ([TJ) was quite arbitrary, and we could have pulled out some 
other factor, say \ s (or A^). This would have affected the Wilson coefficients 
at the matching scale, but the final result would have remained, as it must, 
unchanged. 

We now compute the branching fraction. We evaluate a s in the MS 
scheme (using a s (Mz) = 0.124 as before) and extend the range down to the 
charm scale using the Bernreuther matching conditions03 at the threshold 
fi = mb = 4.87 GeV. Note that we have also taken the CKM matrix elements 
to be real and have neglected any possible imaginary components. Given the 
small values of Ct^{Mw)i this approximation is well justified. It is clear from 
the above that the C2(Mjy) term completely dominates, due to the small 
contributions to 07^ (Mw) from the light internal quark masses. This is in 
stark contrast to b —>■ s transitions (and likewise to s — > d), where the heavy 
internal t-quark forces the the magnetic dipole coefficients to be competitive 
with C2{My/). This can be seen explicitly by comparing the above with 
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Eq. flTEp. The QCD-corrected branching fraction is then 



B 



6a _ V; b V ub 2 

7T V cs 



g(m s /m c ) + 




^g{m d /m c ) 



V 1 



Vj 1 



(4.21 - 7.94) x 10 



-12 



(26) 



where the lower (upper) value in the numerical range corresponds to the 
scale fj, = 2m c (m c ). We see that the effects of the QCD corrections are quite 
dramatic in charm radiative decays, and that the rate is given almost entirely 
as a consequence of operator mixing. The stability of this result can be tested 
once the complete next-to-leading order corrections to the magnetic dipole 
transitions are known. 

Finally, we wish to comment further on the CKM dependence of the 
|Afe| = 1 and |Ac| = 1 effective hamiltonians. We consider each case sepa- 
rately: 

(i) \Ab\ = 1 transition: Here, the the t-quark contribution is seen to 
dominate in every respect. Thus, for the dipole operators Oj and 0$, the u- 
quark and c-quark loops are omitted because they are numerically tiny (e.g., 
see Table 2). Likewise, due to the smallness of the u-quark CKM factors, 
the approximation is made in the literature© to omit any current-current 
operators containing w-quark fields. This explains why only the c-quark 
dependent operators 0\ 2 appear in the |A6| = 1 operator basis of Eq. fllpp . 
This assumption also explains another aspect of the analysis. Ordinarily one 
would expect 0\ 2 to be accompanied by the CKM factor — A c , yet it is the 
prefactor X t which appears in the effective hamiltonian of Eq. @. This is 
because the tiny value of X u has allowed one to write the CKM unitarity 
relation as A c ~ — \ t and thus remove dependence upon X u . 

(ii) |Ac| = 1 transition: In this case, the CKM dependence is more 
complicated since no single quark- loop is dominant. One must expand the 
operator basis as we did in Eq. ([[£]). However, we wish to take note of a 
seemingly remarkable feature which occurs upon carrying out the RG analy- 
sis. The operators 2a and 2 b turn out to have equal anomalous dimensions 
and thus c 2a and C2& have the same numerical coefficient in Eq. (^). The 
most elegant way to understand this result is to exploit the [/-spin symmetry 
present in the system of operators Oi a ib and 2a ,2b- Thus, suppose instead 
of proceeding as we did, we replaced the operators 2 a and 2 b of Eq. (|I~9]) 
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with the equivalent set 0' 2a and 0' 2b , where 

2a ^O' 2a = \{0 2a + 2b ) 

= ^ {u a l^PLS a ){s^P L cp) + {u a ^^P L d a )(d^P L c p ) 

2b ^O' 2b = ^(0 2a -0 2b ) (27) 

= - {.u a ~i^P L s a ){s^P L cp) - (u a jpP L d a )(dp7i*P L Ci3) 

along with a corresponding replacement of coefficients, 

c 2a -> c 2a and c 26 -> c' 26 . (28) 

The matching conditions for the modified coefficients would then have be- 
come 

4,(^ = ^-^ = 1 and c' 2b (M w ) = ^-^ , (29) 

— Aft —A;, 

with analogous replacements made also for Oi a and Oib. Since the mixing 
of 2a and 2 f> with 7 has no dependence on the mass of the internal 
s and d quarks, the operator 0' 2b does not contribute to the process c — > 
m + 7 due to cancellation between the s-quark and <i-quark contributions. 
This cancellation is in fact just the manifestation of an underlying [/-spin 
symmetry. That is, in the limit of neglecting quark mass, 0' 2b carries [/-spin 
1 and thus cannot couple to a photon. This decoupling occurs via the very 
s-quark and <i-quark cancellation under discussion. As a consequence, the 
other operator 0' 2a must have the the same anomalous dimension as 2a 
under RG flow, and we obtain the result cited above. 

As a corollary, it is clear that in the |A6| = 1 transition, the approxima- 
tion made of omitting the w-quark current-current operators is quite unnec- 
essary. One could just as easily deal with an expanded operator basis con- 
taining w-quark fields analogous to that of Eq. (|l9| ) or by invoking an SU (4) 
version of w-quark/c-quark [/-spin symmetry and proceeding as above. 



3 Long Distance Pole Contributions 

As seen in the previous section, the short distance contributions to charm 
radiative decays give very small branching ratios, even when the large QCD 
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enhancement is taken into account. There will, however, be additional con- 
tributions that appear as exclusive modes for which the momentum scale of 
the intermediate quarks is a strong interaction scale and not the short dis- 
tance scale M w . This forces us to view the intermediate states as hadrons 
rather than quarks. These long distance contributions can be partitioned 
into two basic classes. The first corresponds, at the quark level, to annihi- 
lation diagrams cq\ — > g 2 g 3 where a photon line is attached to any of the 
four quark lines. In terms of hadronic degrees of freedom, these give rise to 
the set of contributions which include the pole diagrams. The second type 
of contribution corresponds to the underlying quark processes c — > qi&Q, 
followed by q~2q — > 7. At the hadronic level, this is the so-called vector meson 
dominance mechanism. We shall discuss the pole amplitudes in this section, 
leaving consideration of the VMD mechanism for Section 4. 

The pole amplitudes are but a subset of an entire class of long-distance 
contributions, including the two-particle intermediate states and proceeding 
to all higher n-particle intermediate states. However, of these the most phe- 
nomenologically accessible are the single-particle or pole terms. The relevant 
diagrams, appearing in Figs. 3(a,b), are seen to fall into either of two basic 
classes. We shall refer to transitions as type I if weak-mixing occurs before 
photon emission, i.e. if the incoming D meson experiences weak-mixing, and 
as type II if photon emission occurs before weak-mixing, i.e. if the final 
state meson is created via weak-mixing. In principle, the intermediate states 
occurring in the type I and type II amplitudes consist respectively of all pos- 
sible virtual spin-zero and spin-one particles. We shall find it practicable, 
however, to take into account only the lightest such virtual particles. 



Figure 3. Pole Contributions 
In analyzing long range effects for flavor-changing D decays, we shall 
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employ the effective weak hamiltonian of Bauer, Stech and Wirbeltii-i (BSW) , 

H w = y= [: ai(ud')(s'c) + a 2 (s' d')(uc) :] , (30) 

v 2 

where the colons denote normal-ordering and d', s' are the CKM-mixed fields 

d' = Vud d + V us s , 

s' = V cs s + Vcd d . (31) 
We shall work in the 2x2 basis of quark flavors, 

v={ V v £"W T222 0975) ■ < 32 ' 

\ Vcd Vcs ) \ — u.zzz u.y/o J 

Specific forms for the Cabibbo-favored and Cabibbo-suppressed hamiltonians 
will be given shortly. The quark fields occur in left-handed combinations, 
denoted by 

(91^2) = ^17^(1 + 75)92 , (33) 

and a 1 ,a 2 are free parameters whose values will generally depend on the 
mass scale being probed. Here, they are determined by fitting to D —>■ 
data®, 

ai{m 2 c ) = 1.2 ±0.1 , a 2 (m 2 c ) = -0.5 ± 0.1 . (34) 

Pole Amplitudes of Type I 

Among the possible exclusive D decays, the most promising for experi- 
mental detection occur in the class of vector meson-photon (V7) final states, 

J D(p)^F(k,A) + 7 (q, ( T) • (35) 

For these, the transition amplitude has the gauge invariant form 

M D ->v 7 = el(k,\)et(q,cT) [A pv {jfp v - g^q ■ p) + iA pc e^ af3 k aPp \ . (36) 

The parity-violating and parity-conserving amplitudes are denoted by A pv 
and A pc respectively, and each carries the dimension of inverse energy. Both 
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amplitudes are generally required because the weak interaction does not re- 
spect parity invariance. The D — > decay rate is given by 

r D ^ = ^ (l^ p T + l-4 pc l 2 ) , (37) 

where q is the decay momentum in the D rest frame, 

|q| = ml ~ ml • 08) 
lHl 2m D 

Which particular combination of the parity-conserving and parity-violating 
amplitudes contributes to the decay process will depend upon the weak- 
mixing amplitude. In principle, a charm meson can mix with a sequence 
of either scalar {S n } or pseudoscalar {P n } mesons. Although some work 
on scalar-mixing has been done, EH the outcome is rather model-dependent 
because detailed experimental and theoretical understanding about scalar 
states is lacking. In this paper, we shall therefore consider only the weak- 
mixing of charm mesons with light pseudoscalar mesons and thus work with 
only parity-conserving (PC) pole amplitudes. 

It is appropriate at this point to comment on the notation to be employed 
from this point on in both Section 3 and Section 4. We shall denote fp as 
the decay constant of pseudoscalar meson P and define hy 1 p as the coupling 
constant for the EM interaction vertex of the photon 7 with the mesons 
V, P. Also, the decay constant of vector meson V is given in terms of the 
V-to- vacuum matrix element of the vector current, 

(0K\V\q, A)) = ^^(q, A) = a-VjCq, A) • (39) 

Note that we define two equivalent parameterizations, gy (with units of 
GeV 2 ) and fy (dimensionless), for the vector decay constant. We have found 
that employing g v in the discussion of pole amplitudes alleviates notational 
confusion which would otherwise occur between the vector and pseudoscalar 
decay constants fy and fp. However, it is traditional to use fy in discussing 
VMD amplitudes, and we do so in Section 4. The constants {fv} are ob- 
tained from Fy^i+i- data and have recently been compiled in Table 1 of 
Ref. Q. 

Now, a pseudoscalar state P which is created by weak mixing will prop- 
agate virtually until it eventually decays into the final state. This latter 
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transition is electromagnetic and hence parity-conserving. It has the ampli- 
tude 

M VlP = h VjP el(k, X)et{q, a)e^k aPf3 . (40) 

The absolute value of the coupling constant hy-yp can be inferred phenomeno- 
logically by using 

2 _ / 127rrW 7 /|q| 3 (M V >M P ) 
lHv " pl ~ \ 47rr P ^ 7 /|q| 3 (M P > My) . (41) 

The general type I decay amplitude Aj for D — > V7 is then given by 

Af{D -V 7 )=£ V 7 p. ■ ~ 2 * , ■ (Pn\H^\D) . (42) 

n m D ~ m P n 

With Fig. 3( dbj clS du guide, the notation should be self-evident. 

Predictions for D — > V^j decay amplitudes will be obtained below in 
terms of both type I and type II pole amplitudes, and in the next section we 
shall do the same by using VMD amplitudes. We can, however, accomplish 
somewhat more. In principle, the discussion for final states extends to 
a larger set of meson-photon final states Mj, where the only restriction on 
meson M is that it have spin greater than zero. For each different type of M7 
final state, there will be a gauge invariant D-decay amplitude like Eq. (|36|) 
and an M7P interaction vertex like Eq. (]40f) . However, the generic form 
of Eq. fl42|) continues to hold, except that hy 1 p n is replaced by hwyp n - Of 
course, to have predictive power requires knowledge of the hMyP n coupling 
constant. Fortunately, much has been learned about radiative decays in light 
meson systems over the years. In particular, there are varying amounts of 
experimental evidence for 17 such transitions in the listing of Ref. 0. Of 
these, 10 involve 1~ — > 0~ mesonic transitions, 3 involve 2 + — ► 0~, 2 involve 
1 + — > 0~ and 2 involve 0~ — > 1~. This information allows us to extend the 
analysis of type I amplitudes from just final states to include both 
and T7 configurations as well, where 'A' and 'T' stand for axialvector and 
tensor mesons respectively. The A7 final states are very analogous to the 
V7 decays in that the coupling constant /iA 7 p n is found via Eq. (|4lD and the 
D — ► A'-f decay amplitude has the same form as Eq. (37). For the T7 final 
states, one uses instead 

^ = 40*^ (43) 
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as well as 

T^t 7 = \AT\ 2 • (44) 

To summarize, we shall include in our study of type I amplitudes certain 
D — > A'-f and D — > T7 transitions. For the type II or VMD amplitudes, 
however, we shall limit our calculations to just the V7 final states. 



Cabibbo- favored (CF) transitions: 

In this case, the BSW hamiltonian becomes 

?4 CF) = ~V ud V cs *^ [: ai(iid)(sc) + a 2 (sd)(uc) :] . (45) 

The calculation of weak-mixing matrix elements of D's with the light pseu- 
doscalar mesons is straightforward and results are tabulated in Table 4. The 
fact that these mixing amplitudes are evaluated in vacuum saturation makes 
the forms in Table 4 easy to interpret. Thus, for example, in Cabibbo-favored 
D+ decay, it is the term in the BSW hamiltonian with coefficient a x which 
contributes, and as such, the weak-mixing matrix element is naturally pro- 
portional to the decay constants f n and f D +. 



Table 4 Cabibbo-favored Mixing Amplitudes 
Mixing Matrix Element 

D o ^ ko a 2 V ud V cs *f K f D m 2 D S G F /V2 

D+ -> 7T+ 



For the decay constants of the light mesons we use 

U = 131 MeV and f K = 161 MeV . (46) 

The present situation for the decay constants fo and fo s of the charm mesons 
is somewhat problematic. Experiment provides the upper limit for f D , 

f D < 290 MeV , (47) 
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as obtained from the branching ratio determination £?rD+_^+„ < 7.2 x 10 4 

at 90% confidence-level®. Thus only theoretical estimates exist for fo- 
These occur in three categories, lattice theoretic©. QCD sum rule© and 
quark model fits to color-hyperfine mass splittings'". Estimates fall in the 
range 185 < f D (MeV) < 262. We shall adopt the value 

yjatt. _ 216 MeV ^ ( 4 g) 



which is an average over the lattice estimates'"' and falls between the other 
two types of determinations. 

Recently, the following experimental results (in units of MeV) for fo s 
were announced by CLEolH, WA75© and BEsil, 

f 344 ±37 ±67 (CLEO) 
f D = I 232 ± 45 ± 52 (WA75) (49) 
[ 434±H| i| (BES) , 

where the CLEO value is inferred from the ratio r^n- +I/ /r_ D +_ >( ^ r+ along 
with the branching ratio B D +^^ n+ . In our numerical analysis, we shall use 
the following weighted average of the above decay constants, 

jgpt. = 29g MeV _ ( 5Q ) 

For the sake of comparison, we note that this value is somewhat larger than 
the central value of a weighted average taken from a compilation of existing 



lattice estimates,^' 



yiatt. = 242 MeV _ ^ 

The only other ingredients needed are the radiative coupling constants tiMyP, 
which were defined earlier. Putting together all the necessary ingredients and 
ranging over the set of final state mesons M = p(770), if* (892), &i(1235), 
ai(1270), a 2 (1320) and K%(143Q) yields the magnitudes of type I pole-model 
amplitudes (in units of GeV -1 ) given in Table 5. 

Table 5 Type I Cabibbo-favored Decay 



Mode 




\A?\ (GeV" 1 ) 


Df^ 


p+7 


8.2 x 10" 8 


Df^ 


67/(1230)7 


7.2 x 10~ 8 


Dt^ 


at (1270)7 


1.2 x 10~ 7 


Dt^ 


4(1320)7 


2.1 x 10~ 7 


D° — > 


i?*°7 


5.6 x 10" 8 
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These values should be considered as upper bounds for the following rea- 
son. We have considered the lightest possible intermediate states, pions and 
kaons, because only for these particles is there sufficient data for determining 
coupling constants. However, the pion and kaon intermediate states propa- 
gate far off-shell. Instead of having a squared momentum near the mass-shell 
value q 2 = m 2 , the virtual pion must carry q 2 = m 2 Dg ^> m 2 and similarly 
for the kaon. This effect could well suppress the transition amplitude. 

In principle, one is to sum over all pion-like and kaon-like intermediate 
states. Other possible contributions should be heavier and thus less affected 
by this suppression effect. For pion-like intermediate states, the next state in 
order of increasing mass would be 7r(1300) and beyond that the unconfirmed 
state 7r(1770). Although there is not sufficient data to make a numerical 
estimate of their effect, we can anticipate for such states that 

(i) the propagator contribution will indeed be larger, 

(ii) the weak-mixing between a ground state D meson and a radially 
excited meson P n will be wave-function suppressed, 

(iii) the radiative coupling constant tiM-yP n might well be relatively smaller 
due to phase space competition with other decay modes of the massive meson 

We would expect the net result of these effects to decrease the overall 
contribution from the excited states. 



Cabibbo-suppressed (CS) transitions: 

The weak-mixing now proceeds according to the weak hamiltonian 



n 



(CS) 



V2 



: ai (V ud V cd *(ud)(dc) + V us V cs *(us)(sc)^ 



+a 2 (V us V cs *(ss)(uc) + V ud V cd \dd)(uc)) : 



(52) 



The action of the (dd) and (ss) operators on the vacuum when expressed in 
terms of the pseudoscalar meson states becomes 



(dd) = -0.70717T + 0.58?7 + 0.40?/ 
(s S ) = -0.5777 + 0.82r/' , 



(53) 
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where an rj — rj' mixing angle Op = —20° is adopted, crl In addition, we 
takeH 

fv - fv> - f* • ( 54 ) 
The mixing amplitudes which are relevant for Cabibbo-suppressed decay ap- 
pear in Table 6. Observe that we have simplified the notation for D° transi- 
tions by expressing V u dV c / in terms of V US V CS * . 



Table 6 Cabibbo-suppressed Mixing Amplitudes 

Matrix Element 



Mixing 


D + -» 






K+ 


D° -> 


7T° 


D° — > 


V 


D° — >■ 





aiV u dV* d f 7T fD+ r nT 2 D+ G F / V2 
axV us V^j K j Ds m 2 D G F /V2 
0.7Q7la 2 V us V*f DO Um 2 DO G F /2 
-1.1 5a 2 V us V*J D o f v m 2 D0 G F /2 
OA2a 2 V us V c *J D of v ,m 2 DO G F /2 



The analysis for Cabibbo-suppressed decays proceeds analogous to that 
for Cabibbo-favored decays, with one significant complication. For each of 
the Cabibbo-favored transitions, only one amplitude contributes. For D° 
decay, however, all the Cabibbo-suppressed pole amplitudes contain a sum 
over 7T°, r) and rf intermediate states. It is important to get the relative 
phases of the interfering amplitudes correct. We have therefore performed 
an analysis of the nine V° — > P°7 couplings in light of the most recent data, 
where V° = p°, to , and P° = ir°, r\ and t/.KU The magnitudes of the 
Cabibbo-suppressed amplitudes are displayed in Table 7. 

Table 7 Type I Cabibbo-suppressed Decay 
Mode pjfj (GeV 1 ) 

1.3 x 1(T 8 
1.2 x 1(T 8 
4.9 x 1(T 9 

3.4 x 1(T 8 
2.8 x 10~ 8 

6.0 x 1(T 8 
4.8 x 10~ 9 

6.1 x 1(T 9 
7.4 x 1(T 9 



D+ - 


-> p+7 


D+ - 


-> 6^(1230)7 


D+ - 


-> a^(1270)7 


D+ - 


-> a+(1320)7 


Dt- 


-> K*+~? 


Df- 


-> Kt( 1430)7 


D° - 


-A 


D°- 




D° - 


+ 0° 7 
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Pole Amplitudes of Type II 

Analogous to the type I D — > decay amplitude of Eq. (K2|) we have 



AfAD - V-y) = £ (V\H^\D* n ) ■ 1 ■ (55) 

for the corresponding type II transition. From the viewpoint of phenomenol- 
ogy, the type II transitions are more problematic than are those of type I 
because less experimental input is available. Thus, we shall need to rely a 
bit more heavily on theoretical predictions. 

The first difficulty is that the couplings h D *o^ D o, /i_d*+ 7 d+, and h D *+ lD + 
have not yet been experimentally measured. This is because, although the 
relevant photonic branching ratios have been measured, only upper bounds 
exist for the full widths of the associated spin-one exited states, D*°, D* + 
and D* + , 

T D *o < 2100 keV =^ r D *o^ D o 7 < 764 keV 
T D *+ < 131 keV =>• T D *+^ D+7 < 1.44 keV (56) 
r D .+ < 4500 keV T D *+_ ¥D +^ < 4500 keV . 

Fortunately, predictions for the T D *o^ D o^, T D *+^ D + y and T D *+^ D + ry transi- 
tions have appeared in the literature recently. H3' There is some 



spread in predictions, and so we choose the representative values, 

r , £>»o_ >jD () 7 = 20 keV , T D *+_> D +^ = 0.5 keV , T D .+^ D + 7 = 0.3 keV , (57) 
which implies 



(58) 



where we have adopted the phases implied by the quark model. A rough check 
on whether the above values are reasonable is afforded by the nonrelativistic 



h D , 


o jD o 


= 0.542 GeV" 1 , 




h 7 £)+ 


= -0.087 GeV 1 


S 




= -0.066 GeV" 1 
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quark model, in which 



h D *o jD o 

h D*+ 1 D+ 



2e 



1 1 

M~ c + W u 



2 
2 



1 

M~ d 
1 



(59) 



where the {M^} are constituent quark masses, distinct from the current 
masses {m^} of Section 2. If we take M c ~ 1.64 GeV, as implied by a fit to D 
and D* masses, then the relations in Eqs. (|57| - |59D yield M u ~ Md — 0.48 GeV 
and M s ~ 0.53 GeV. 

The other of the difficulties concerns the weak-mixing matrix elements. 
For type II transition amplitudes, the mixing occurs between charm and light 
vector mesons, as in 



(60) 



In the above, the gy are the vector meson decay constants defined in Eq. ([39]) 
and whose determination we shall discuss shortly. As with the type I am- 
plitudes, we have employed vacuum saturation. To determine the action of 
the (dd) and (ss) operators upon the vacuum we employ the ideally-mixed 
vector meson states, so that 



(dd) 



u - p u 



and 



{ss 



(61) 



For the light 1~ mesons, the collection {gy} of vector decay constants can 
be determined by referring to the vacuum-to-meson matrix elements of 
given in Table 1 of Ref. [[|. Together with isospin and SU(3) relations along 
with quark model insights, these generate all the needed values, e.g. 



9 P ^ 



0.17 GeV 2 



gK* 



2 



m. 



0.22 GeV 



(62) 



To estimate the D* + and D*° decay constants, we invoke the heavy-quark- 
symmetry relations, 

0.588 GeV 2 , 



9d* 
9d* 



m D f D ~ 0.403 GeV 2 



(63) 
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The magnitudes of the type II amplitudes thus calculated are given in Table 8. 



Table 8 Type II Decays 



Mode 




l^zl (GeV-) 




■ p+7 


1.9 x 10" 8 


D° — > 




5.9 x 1(T 8 


D+ -» 


' P + 7 


3.6 x 1(T 9 




. k*+ 7 


5.1 x 1(T 9 


D° — ► 


P°7 


4.7 x 1(T 9 


D° — ► 




6.9 x 1(T 9 


D° — ► 


</>°7 


1.6 x 1(T 8 



D* excitations with spins not equal to one will not contribute to type II 
amplitudes if we adhere strictly to the hamiltonian of Eq. (pO|)^l and continue 
to work within the vacuum saturation framework. The reason is that mesons 
with J > 1 cannot have a nonzero matrix element with the vacuum via the 
current 97^(1 + 7s)c. The possibility of an intermediate charm meson with 
J = is disallowed since it could only mix with a final-state J = particle 
and the decay of a spinless particle to another spinless particle plus a photon 
is forbidden. 

Although we have considered just final state vector mesons in Table 8, 
it should be obvious that in principle the spin-one intermediate D* states 
can also mix weakly with axialvector final state mesons. Unfortunately, one 
knows less about the decay constants of axialvector mesons than one does of 
the vector mesons. 

4 Long Distance VMD Contributions 

The VMD contribution to charm meson radiative decay is depicted in Fig. 4, 
where a D meson is seen to (i) decay weakly into a final state of a vector meson 
V and a meson M of nonzero spin, followed by (ii) an electromagnetic VMD 
conversion of V into a photon. Roughly speaking, in the VMD approach the 
D — > M7 amplitude is obtained by multiplying the D — ► MV amplitude 
by the factor e/ fy where e is the electric charge and fy is the dimensionless 



24 



version of the vector meson decay constant defined in Eq. (|39|). It is important 
to keep in mind that in the VMD process D — > MV, the vector meson V is 
off-shell. Thus, to obtain the VMD amplitude for D — > M7 will require an 
extrapolation from = to p y = for both the V — > 7 vertex and the 
D — > transition. For our considerations, the main intermediate states 
will involve virtual rho and phi mesons. We shall employ the observation 
made in Ref. that the rho-gamma vertex seems to be unaffected by 
the extrapolation whereas the phi-gamma vertex is reduced by a factor of 
77^ ~ \J~2. In the following, we will consider a number of examples for the 
case M = V, and so we shall be working with VMD chains which begin with 
the process D — > VV. Since the VV final states have L = 0, 1, 2 as allowed 
orbital angular momentum values, the VMD amplitude will in general have a 
parity-conserving part ^4 V md corresponding to the VV P-wave and a parity- 
violating part -4-vmd corresponding to the VV S- wave and/or D-wave. 



Figure 4. VMD Contribution 

In practice, there are two means for determining the D — > MV part of a 
VMD amplitude for D — > M7: 

1. One can input D — ► MV" experimental data directly in order to phe- 
nomenologically determine the D — ► MV" amplitude. In this approach, 
it is crucial to maintain gauge invariance. A careful discussion of 
how to construct a gauge invariant amplitude was recently given in 
Ref. (which considered this type of empirical VMD contribution 
to B — > K*j), so we need not detail this procedure here. Since the 
database for D — > MV transitions is unfortunately small, the ability 
to generate VMD amplitudes using this phenomenological method is 
limited. 
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2. One can employ some theoretical description to model the D — > MV 
amplitude. Since the models currently available do not always reli- 
ably reproduce branching ratios and polarizations of final-state vec- 
tor mesons in decays of heavy mesons,® this method is also not 
beyond criticism. For definiteness, we shall continue to employ the 
BSW modelEB introduced in Section 3. Within this approach, the 
squared VMD amplitude for the important case where M is a vector 
meson becomes 



VMD 



G 2 F\V*qVqu\ 

2m 2 D k 2 



(m D + m Y y A{(q 2 ) + 



Ak 2 m 2 D V 2 (q 2 ) 



x 



v fv , 



(64) 



where |k| is the photon spatial momentum, q represents either of the d 
or s light quarks, and J is a process-dependent isospin coefficient. The 
BSW coefficients ai{m 2 ) and a 2 (ml) which correspond to the color- 
favored and color-suppressed operators are given in Eq. (|34"1) . The 
remaining notation is explained by noting that in the factorization 
approximation for D — ► MV, one of the final state particles, which 
we call X (either M or V), couples directly to the vacuum and the 
other, which we call Y (either V or M) , appears in the D-to-Y matrix 
element of the charged weak current J^ h . Thus the quantity fx is the 
decay constant of X, and Ai(q 2 ) and V(q 2 ) are the semileptonic form 
factors defined by 

(Y( P y)\JiL\D(P)} = 
2V{q 2 



mr> + my 
+i 



e^elP pPYa + 2m Y tA (q 2 ) e -^q» 

Q 

m D + m y )A 1 (g 2 )e^ - ^^1(P + Vy f 



mr> + vti Y 



-2m Y A z {q 2 f-^q' 
q z 



(65) 



In the VMD amplitude, the form factors are to be evaluated at q 2 = 
if X = V and at q 2 = m 2 M if X = M. Throughout, we shall make 
use of the form factors as measure dm in D — > K*lv and also employ 
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SU(3) relations as needed. This should provide a good estimate of 
the form factors appearing in the D-to-p and D-to-0 matrix elements. 
Whenever the form factors are to be evaluated at momentum transfers 
other than at q 2 = 0, we shall use a monopole form to extrapolate from 
q 2 = 0. This amounts simply to dividing the form factors at q 2 = by 
the quantity 1 — q 2 /m 2 ole . 

In the following, we shall give VMD predictions for a number of specific 
D — > M7 decays, grouped as Cabibbo favored, singly suppressed or doubly 
suppressed. In the few cases where we can employ both the above approaches, 
we shall refer to them respectively as 'Meth. 1' and 'Meth. 2'. Given the 
lack of abundant D — > MV data, however, we shall be forced to adopt the 
theoretical approach of Meth. 2 in most cases. 

Before we can proceed, there is another topic which must be addressed, 
the dynamical complication of significant Final State Interactions (FSI). Al- 
though presumably not a problem in B — > MV decays, detectable FSI are 
known to exist in the D-meson mass region. This can produce an ambiguity 
in the VMD analysis because FSI will inherently be part of any VMD am- 
plitude obtained from D — > MV data, but will not be present in the BSW 
construction. It is difficult to remove the effect of FSI from the phenomeno- 
logical VMD amplitude because the vector meson V is to be taken off-shell, 
and FSI might have an important kinematic dependence, i.e. the py depen- 
dence of the FSI has also to be taken into account. Consequently, any FSI 
effects entering in data may not be present to the same extent in the VMD 
amplitudes. As regards the factorization construction (Meth. 2 above), the 
exclusion of any FSI effects in the BSW amplitude amounts de facto to a 
specific prescription for the py-dependence of the FSI. There is some infor- 
mation on the py-dependence of the 7 — V couplings and of certain matrix 
elements, but it is not possible at this time to separate the two effects. As 
we show in the Appendix, the effect for p emission in A 2 (1320) decay can 
be as much as a factor of two. By contrast, no such suppression is seen in 
p-photoproduction, although in 0-photoproduction an effective reduction of 
about y/2 in amplitude is observed and a somewhat smaller effect of vTo~ is 
seen in cu-photoproduction. 

Cabibbo-favored Modes 
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D° — > K°*7: This is an instance in which the phenomenological approach 
is applicable since experimental information on the D — > MV intermediate 
state is available. There is a branching ratio determination© B D o_+K*o p o = 
(1.6 ± 0.4)% and the amplitude is known to be (i) almost all transverse and 
(ii) almost all S'-wave. [] This allows us to write the VMD contribution to 
D° -> iT*°7 as 

4 ¥ md = t ■ aM ; 0p ° , A%»^0, (Meth. 1) (66) 

where we follow the notation of Ref. || and denote the phe- 

nomenological S-wave amplitude for D° — ■> K*°p°. With T d o^k*°p° = 

2.53 x 

1(T 14 GeV and a D o^ R ,o p o = 1.63 x 1(T 6 GeV, thisyields A^ MD (D° -> K*°j) 
of about 6.8 x 10~ 8 GeV -1 . The data on D° — > K*°p° is consistent with no 
parity-conserving (P-wave) contribution. 

Alternatively, the factorization approach of Eq. ( |64l) predicts both am- 
plitudes. In this case, we take aj = a,2, and the vector meson to be mixed 
with the photon is the p°, so that V = p° and X = K°*. The form factors 
needed are those entering in D — > p semileptonic transitions. Making use 
of the measured D — > K* form factors implies 7=1/2. To extrapolate 
the form factors from q 2 = to q 2 = m 2 K *, we use a monopole form where 
the D* is the nearest singularity. The parity-violating and parity-conserving 
amplitudes are given in Eq. (Q) by the terms involving the A\ and V form 
factors respectively. Using fx* = 0.2 GeV 2 we obtain 

^p v md = 5.1 x 10~ 8 GeV -1 , ^ C MD = 3.8 x 10~ 8 GeV^ 1 . (Meth. 2) 

(67) 

We notice that -4vmd i s i n reasonable agreement with the one obtained from 
the use of data from the nonleptonic mode, given the large uncertainties 
involved in these predictions. Indeed, the factorization estimate for the D° — > 
K°*p° S'-wave amplitude gives a D o^x*°p° = 1-3 x 10~ 6 GeV -1 which is within 
20% of the experimental value. It also predicts a P-wave branching fraction 
of 0.15% for D° — > K°*p°, which is below the current upper limit of 0.30%. 

1 The Particle Data Group also lists branching ratios of (3.0 ± 0.6)% and (2.1 ± 0.6)% 
for S'-wave and D-wave respectively. H These values are completely consistent with the 
fact that the total transverse mode (which must be entirely S'-wave by the absence of any 
P-wave) is (1.6 ± 0.5)% and that the S-wave (longitudinal) must cancel with the D-wave 
to produce the net zero longitudinal branching ratio. 
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— > p + 7^ The VMD amplitude for this decay proceeds via Df — > 0p + 
followed by 0-7 mixing. Although the branching ratio for Df — > 0p + is 
known to be (6.5 ±J"| %), no information on helicities or partial waves exists, 
so we cannot apply the phenomenological method here. Turning instead to 
the factorization approach of Eq. (|64D , we have X = V = and Y = M = p + . 
Therefore we require the — > <fi semileptonic form factors evaluated at 
<?o = m p- Although there is experimental information on these decays, the 
branching fraction and the form factors depend strongly on -B D +_^ 7r +, which 
is still very uncertain. Thus, again making use of D — > K* data, taking 1 = 1 
and with a decay constant of g p ~ 0.17 GeV 2 , we find 

^p v md = 3.2 x 1(T 8 GeV" 1 and ^ C MD = 2.8 x 1(T 8 GeV" 1 . (68) 

The Cabibbo-favored VMD amplitudes are summarized in Table 9. 



Table 9 Cabibbo-favored VMD Amplitudes 



Mode 


^VMD (10' 


- 8 GeV" 1 ) 




Parity- conserving 


Parity- violating 


D° — > 




3.8 


5.1-6.8 




P + l 


3.2 


2.8 



Singly Cabibbo-suppressed Modes 

D° —>■ p°7: This process can proceed via two different intermediate states, 
namely p°0 and p°p°. There is one known branching ratio B D o^ p o,p = 
(1.9 ± 0.5) x 10~ 3 with no helicity (or partial wave) information. Letting 
t)t be the transverse fraction of the observed branching ratio, rjs the 5-wave 
fraction in the transverse mode, and r\p the P-wave fraction, we then obtain 
for the S"-wave amplitude of D° — > p°0, 



a D ^ p o, = m D ^ TD0 ^ Tr]S = 7 x lO" 7 ^ GeV , (69) 
and for the corresponding P-wave, 



m^rrip y |k| 3 



4^0 o, VTVP = iM x 1Q - 6v ^ GeV -1 , (70) 
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where again we employ the notation of Ref. in denoting b D o_^ p o^, as the 
phenomenological P-wave amplitude for D° — > p°<f). Then, multiplying by 
the VMD factor e/f^, we obtain the Meth. 1 estimate 

^p v md = 0.60 x 10" 8 GeV" 1 and *4^ MD = 1.0 x 10" 8 GeV (71) 

for i]t ~ 0.5 and r]$ ~ 0.66. On the other hand, there is no available 
experimental information for D° — > p°p°, other than B D o_^ 7r + n - 7r + n - = (8.3 ± 
0.9) x 10~ 3 which can be taken as an upper limit. Let us also estimate the 
D° — > p°7 mode in the factorization approach, which can be used to predict 
both the off-shell amplitudes, D° — > p°<f) and D° — > p°p°. In both cases we 
need the D° — > p° form factors, for which I = 1/2. Using Eq. (^) we obtain 

^vmd(pV^A) = 0.22 x 10- 8 Gev" 1 

^vmd(pV P°7) = 0.18 x 10- 8 Gev" 1 

^vmd(pV ^ P°7) = 0.75 x 10- 8 Gev" 1 

^vmd(pV-P°7) ^ . (72) 

Our estimate for the parity-conserving p°p° — > p°7 transition is based on the 
observation that an on-shell P-wave p°p° state is forbidden by Bose statistics 
and hence the associated off-shell amplitude will be suppressed. First, let us 
compare the first two rows in Eq. (|72]) with the results obtained in Eq. ([H]) 
by making use of the p°0 data. We can see that the factorization amplitudes 
are lower, as caused by smaller predictions for the nonleptonic intermediate 
modes. In general, factorization predictions will be modified by FSI. For 
instance, in the case at hand there could be a large enhancement due to 
K*K* — > p°(fi rescattering effects, t^l If this is the case, this effect strongly 
depends on the kinematics and it is different in the off-shell nonleptonic am- 
plitudes entering in the calculation of the VMD diagrams. The factorization 
approach provides a prediction which is free from FSI effects. In these cases 
we will take these two estimates as the allowed range. On the other hand, 
factorization predicts that the p°p° intermediate state provides most of the 
VMD amplitude. When both intermediate states are taken into account in 
the factorization estimate, the predictions of Eq. (|7lD and Eq. ( |72|) roughly 
agree. This will not be the case for the following mode. 
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D° — > 07: Now, there is only one nonleptonic intermediate state, 
The amplitudes as extracted from D° — > <pp° data are 




On the other hand, factorization predicts the much smaller amplitudes 



^p v md = 0.7 x 1(T 8 GeV" 1 and ^ C MD = 0.6 x 10" 8 GeV" 1 . (74) 



Part of the difference between the predictions in Eq. (|73|) and Eq. ([74]) may 
be due to the presence of FSI effects in the on-shell amplitude measured and 
used in Eq. fl73]) and the assumed absence of FSI in Eq. (f74|) . In Table 10 we 
include both predictions as the allowed range. 

D° — > w-f: This mode is very similar to D° — > p°7 and we obtain from the 
factorization approach of Meth. 2 



^p v md = 0.7 x 10" 8 GeV- 1 and ^ C MD = 0.6 x 10~ 8 GeV" 1 . (75) 



^p c md = 1.9 x 10~ 8 GeV" 1 and ^ V MD = 1.6 x 10~ 8 GeV" 1 . (76) 

Incidentally, the expectation for D + — > p + p° is that its branching ratio should 
be at least 0.4%. 

Df — > K* +/ ~f. Proceding analogously, we use the factorization estimate of 
£>+ — ► K* + p° to express the VMD amplitudes for Df — > K* + >y as 



D + — > p +r y: Here, the mode D + — > p + p° should give the dominant contri- 
bution to the VMD amplitude, implying the Meth. 2 amplitudes 



-VMD 




(m Ds + m K *) 



At 



m Ds E, 



'7 ' 



(77) 



A P c 



-VMD 




-VMD ' 



2 V 

(m Ds + m K *) 2 A 1 



m K *m p 
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Upon taking the form factors to be identical to those in D + — > K*° in the 
SU(3) limit, we have 

^p v md = 1.0 x 10~ 8 GeV" 1 and ^ C MD = 0.9 x 10~ 8 GeV" 1 . (78) 

Our VMD predictions for the magnitudes of the Cabibbo-suppressed tran- 
sition amplitudes are summarized in Table 10. 



Table 10 Cabibbo-suppressed VMD Amplitudes 



Mode 


I^vmdI (10- 


8 GeV- 1 ) 




Parity-conserving 


Parity- violating 


D+ -» 


p+7 


1.6 


1.9 


Df^ 




0.9 


1.0 


D° — > 


A 


0.2 - 1.0 


0.5 - 1.0 


D° — > 




0.6 


0.7 


D° — > 


<A 


0.6-3.5 


0.9-2.1 



Doubly Cabibbo-Suppressed Modes 



Finally, to estimate the size of the doubly Cabibbo-suppressed modes, we 
consider the D + — > K* +, ~f transition. Upon computing the amplitudes using 
the factorization expression of Eq. (|64j) , we obtain 

•^vmd = 4.2 x 10~ 9 GeV" 1 and ^ C MD = 4.4 x 10~ 9 GeV" 1 . (79) 

Similarly, we find for the mode D° — > K* 0, y 

A^ MD = 1.75 x 10~ 9 GeV" 1 and ^ C MD = 1.83 x 10~ 9 GeV" 1 . (80) 

5 Summary and Conclusions 

The existing database for direct evidence of radiative D decays is meagre, 
as has been shown in Table 1. However, interesting levels of experimental 
sensitivity are currently being attained and we can anticipate the detection of 
radiative signals in the not-too-distant future. Our motivation in undertaking 
the study reported here has been to stimulate such experimental efforts. 
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As shown here, weak radiative decays of charmed mesons are not domi- 
nated by the short distance penguin diagrams of Figs. 1(a), 1(c), but rather 
by long distance processes involving nonperturbative strong interaction dy- 
namics. From the standpoint of probing the inner workings of the Standard 
Model, one might have naively hoped to use radiative charm decays in order 
to observe the short distance process c — > wy. This would be in analogy 
with radiative .B-meson decay, where the amplitude is dominated by the 
penguin transition b — > s'y and receives a large enhancement from QCD ra- 
diative corrections. By contrast, the corresponding c — > wy charm transition 
is minuscule at lowest order and would require an unexpectedly large QCD 
enhancement to become detectable. To our knowledge, the calculation per- 
formed here of the QCD radiative correction to c — > wy is the first explicit 
and detailed analysis of this system given in the literature. In addition, we 
were able to employ [/-spin arguments to clarify the role played by neutral, 
flavor- changing operators such as ucqq which contribute to the expanded 
operator basis in the RG analysis. Our conclusion that the c — > wy QCD 
radiative corrections are substantially larger than for b — > s'y is due in part to 
the large operator mixing at the lower renormalization scale associated with 
the c-quark and in part to the disparate sizes of the Wilson coefficients at the 
matching scale of the contributing operators. Nevertheless, the radiatively- 
corrected c — > wy penguin transition remains extremely small. The main 
sources of suppression are the small quark masses and also the CKM factors 
iKfe^l 2 occurring in the numerator of the c — > wy branching fraction in 
Eq. fl26|). In the Wolfenstein parameterization, with sin# c = A, this CKM 
dependence amounts to a A 10 suppression in decay rate. 

On the other hand, we have shown in Sections 3 and 4 that long distance 
contributions are several orders of magnitude larger. A very rough estimate 
of the typical branching ratio to be expected is 

Bd->Mj ~ at cm B™^M Pt ' , (81) 

where B^^J pi " is the branching ratio for the nonleptonic D decay to some 
final state M. Thus, typical branching ratios of order B^^J pt " ~ (0.001 — 
0.05) would induce radiative branching ratios in the range Bn^M-y ~ (7 x 
10~ 6 — > 4 x 10~ 4 ). In Sections 3 and 4 we have performed a more detailed 
analysis by modeling the nonperturbative dynamics. 

Inspecting the long distance contributions to the set of exclusive processes 
D° — ► p7, D° — > wy, D s — > K* +, y and D + — > p +/ y for which c — > wy is 
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the underlying transition, we see that the VMD and pole amplitudes carry a 
single factor A, therefore representing an enhancement of A 4 over the penguin 
amplitude. As can be seen in Table 11 the expected branching fractions for 
these modes are in the 10~ 6 — 10~ 4 range, whereas we estimate B C ^ U1 ~ 10~ 12 . 
As a consequence, c — > wy is not a good process to test the validity of the 
Standard Model. That is, a hypothetical contribution from new physics 
would have to be extremely large in order to overcome the long distance 
physics. 

The situation is very different in radiative B decays. The short distance 
transition b — *> S7 has the same CKM structure as the corresponding long 
distance contributions. For instance, the mode B — > K*y might conceivably 
have long distance contamination of the order of 20% in the rate 0. Although 
this is small compared to the charm case, it would be desirable to reduce the 
uncertainty in the calculation of these effects in order to subtract them from 
the measured signal. Moreover, long distance effects could also be affecting 
the inclusive b — > 57 branching ratio, therefore limiting the precision with 
which the Standard Model can be tested in these decays. 

The various amplitudes are summarized in Table 11 and are given there in 
units of 10 -8 GeV~ . In principle, the most conservative attitude is to take all 
relative signs as unknown, which would render the calculation of branching 
ratios highly uncertain. Fortunately, with the aid of the quark model we can 
reduce this ambiguity. The relative sign of Pole-II to Pole-I contributions is 
affected by (i) a minus sign difference in the pole denominators, (ii) an extra 
minus sign in type-II amplitudes due to the vector meson propagator, and 
(hi) minus sign differences in the VPy couplings between the c-quark and 
light-quark EM sectors. In a quark description of a qiq2 meson, this latter 
sign is inferred by studying 



Although this line of reasoning narrows down the range of predictions sig- 
nificantly, experimental data will be needed to obtain information regarding 
the relative phase between the pole and VMD contributions. In this regard, 
it will be helpful to note that, at least in our approach, the parity- violating 
amplitudes arise solely from the VMD process. 




(82) 
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Table 11 Amplitude and Branching Fraction Predictions 



Mode 






A pc 


A pv 


B njM , flO" 5 ) 




P-T 


P-TT 

J. .L.L 


VMD 


VMD 




D + n + 'y 


8 2 


— 1 9 


±3.2 


±2.8 


6 — 38 


D°K*°-f 


5.6 


-5.9 


±3.8 


±(5.1 - 6.8) 


7-12 




7 2 








~ 6 3 




1.2 








~ 0.2 




2.1 








~ 01 




1.3 


-0.4 


±1.6 


±1.9 


2-6 




1.2 








~ 3.5 




0.5 








~ 0.04 




3.4 








~ 0.03 


D+K*+7 


2.8 


-0.5 


±0.9 


±1.0 


0.8-3 


DtK* 2 + 1 


6.0 








~ 0.2 


D°A 


0.5 


-0.5 


±(0.2-1.0) 


±(0.6 - 1.0) 


0.1-0.5 


D°cA 


0.6 


-0.7 


±0.6 


±0.7 


~ 0.2 


£>V°7 


0.7 


-1.6 


±(0.6-3.5) 


±(0.9-2.1) 


0.1-3.4 




0.4 


-0.1 


±0.4 


±0.4 


0.1-0.3 


D°K*°-y 


0.2 


-0.3 


±0.2 


±0.2 


~ 0.01 



Finally, let us comment on the inclusive photon spectrum. In the B 
system, the quark transition b — > S7 provides a useful framework for pre- 
dicting properties of the hadronic inclusivedecay B — > X s 7. Thus, one 
estimates the S — > X s 7 decay rate by computing the b — > 57 decay rate 
and normalizing relative to the semileptonic decays to eliminate undue de- 
pendence on the mass m&. Likewise, one predicts the photon energy spec- 
trum in B — > X s ~f decay by referring to the underlying two-body b — > 57 
decay. H3'0 If quarks were free, there would be a monochromatic photon 
spike at E 1 = (ml — m 2 s )/2mb ~ va^jl. In reality, the photon spectrum be- 
comes broadened via hadronization of the s-quark jet. The individual strange 
mesons (if* (892), if 1 (1270), etc) which populate the inclusive final state X s 
originate predominantly from the s-quark jet hadronization. These explana- 
tions of B — > X s 7 inclusive decay are in accordance with the spectator model, 
and so isospin symmetry should manifest itself event-by-event. For example, 
the rates for isospin-related modes such as B° — > K *°7 and B~~ — > if*~7 
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should be equal. A deviation from this pattern would constitute evidence for 
either non-spectator or new-physics contributions. In a heavy-quark effec- 
tive theory description, such non-spectator effects would occur at subleading 
level. 

The theoretical description of charm inclusive decay could hardly be more 
different. Now, there is no emergent light-quark jet which hadronizes to form 
the set of final states. Instead, the 'black box' of long-range effects such 
as pole-amplitudes, VMD-amplitudes, etc dominates the physics. Thus, to 
determine the photon energy spectrum in D — > X„7, one would sum over 
the most important of the exclusive radiative modes. Presumably this would 
yield a reasonable description at least over the part of phase space where the 
photon energy is largest. It would be prudent to be on the lookout for the 
unexpected. For example, exclusive modes in light meson radiative decay are 
known to exhibit rather large isospin-violating effects, as in 

— ^— — ^— ^ = 2.27 ± 0.30 and I>z^X = 1.76 ± 0.49 . (83) 

If this effect were to be maintained mode by mode in the exclusive D decays, 
it would lead to interesting levels of isospin violation in the inclusive decay. 
Of our results, Table If indicates that the likeliest possibility for isospin 
violation would appear to be in D° — > p 0r y/D + — > p +r y. 

Charm radiative decays give us the opportunity to study various aspects 
of long distance dynamics. We have seen that the theoretical predictions of 
the branching fractions are, in some cases, rather uncertain due to model 
dependence. Experimental information will therefore be needed to complete 
the theoretical picture of these decays. It is an interesting irony that the 
understanding gained from future observation of different D radiative decays 
can then be used to predict more confidently the size of such effects in B 
decays. 
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Appendix: VMD in Light-Meson Radiative Decay 

The original application of VMD for analyzing hadronic radiative decays 



occurred in the light meson sector, e^) In order to test the VMD method using 
an up-to-date database, we too shall consider (briefly) light meson radiative 
decays in this Appendix. As we shall see from our study of two particularly 
clean examples, the situation is encouraging but not uniformly so. First, we 
shall revisit the original arena for testing VMD, the p and uj decays into pion- 
photon final states. Then we shall analyze decays of a higher mass state, the 
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tensor meson ^4^(1320). We stress that in each of these cases the transition 
is purely electromagnetic, unlike the more complicated electroweak decays 
treated in the main body of the paper. Therefore, the 'pole' amplitudes do 
not occur here since there is no weak mixing, so one obtains a clean look at 
the VMD contribution. 

Radiative Decays of the Vector Mesons p and u 



There are three electromagnetic P-wave decays in the p - u system, 



CO — > 7T 7 



P 



7r + 7 and p° 



(84) 



In the VMD approach, these are described in terms of two electromagnetic 
mixing amplitudes, cj-7 and p-7, and one strong interaction vertex, g Wfm . 

Due to the off-shell nature of the VMD amplitudes, different momentum 
regions occur in the oupir vertex for the transitions of Eq. ([84]) . In to — > 717, the 
intermediate p propagates at q 2 = whereas for p — > 777 it is the intermediate 
lo which propagates at q 2 = 0. Part of the VMD folklore built up over the 
years is that extrapolation of the light vector meson squared-momenta from 
the meson mass-shell to the photon mass-shell does not strongly affect the 
decay amplitude. The ratio of p and u decay widths can be used to test this 
as follows. Recall that for the VMD description of 1~ — ► 0~7 transitions, the 
strong vertex is related to the decay width F via 



9i 



127rr 

|q| 



nl/2 



(85) 



where / — > f p for to decay and / — > f u for p decay. Noting that the decay 
momenta in ui — > n'-f and p —> rcy are almost equal, one has 



v, 



3.24 ±0.19 



(86) 



->7T" r 7 



provided the same strong vertex is used in each decay. In the above, we have 
used the charged-p decay width in view of its superior accuracy. The value 
appearing in Eq. (|8E|) is seen to be in accord with that inferred from vector 
meson decay into lepton pairs ( cf Table 1 of Ref . || ) , 



fp 



3.39 ±0.10 



(87) 
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Alternatively, one can use each of these radiative decays to extract deter- 
minations of g ufm as in Eq. (jS5f), and one finds 

f (11.73 ± 0.35) GeV- 1 (to -> 7r° 7 ) 
g uim = \ (12.40 ± 0.64) GeV- 1 (p+ - 7r+ 7 ) (88) 
[ (16.40 ±2.1) GeV" 1 (p° -> 7r° 7 ) . 

The a; — > 7r°7 and p + — > 7r + 7 determinations are seen to be consistent within 
experimental error. This is significant because these decays involve different 
momentum extrapolations as discussed above. The larger coupling obtained 
from p° — > 7T°7 decay has substantially larger errors. We now turn to a 
different transition in which, if one accepts the data at face value, a non- 
negligible momentum dependence is present. 

Decays of the Tensor Meson A£(1320) 

The meson A^"(1320) has been observed to decay into both the irp and 7r 7 
modes, with branching ratios 

&M-+kp = °- 701 ± °- 027 and B^a-wr = (2.8 ± 0.6) • 10" 3 . (89) 

These data turn out to provide a particularly clean test of the VMD method 
in two respects. First, there is just a single partial wave in the final state. As a 
consequence, the decay rates alone can be used to test VMD without any need 
for polarization information of the final state particles. The occurrence of a 
single orbital angular momentum in the final state follows from conservation 
of parity and of angular momentum. Thus we have 

V: + = R 2 R L L = 0,2,4,... (90) 

J : |2| = |1 + L| L = 1,2,3 (91) 

which implies that L = 2. In addition, of the three light vector mesons p,cu,(f), 
only the p can appear together with a pion in a final state of A 2 decay. The 
reason is that the decay A 2 — ► rr V (V is a vector meson) proceeds through 
the strong interactions and conservation of G-parity forbids the ttuj and 7r0 
modes. Thus, the rho is the only light vector meson involved in the VMD 
determination and interference with uj or <p mediated processes is absent. 
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The amplitude for the transition A^ip) — > 7r(q) + p (k) can be written as 

A P = ^ ua %el(k)q a hp a (p)q' J , (92) 
m A 

where g np is a dimensionless quantity and hp a (p) is the spin-two polarization 
tensor of the A?,- From the decay rate relation, 

one determines a magnitude for the coupling g np . This can be used, in turn, 
to predict the radiative coupling g ni via the VMD formula 

<7™ - ~7~9"Rp , (94) 

and we find 

g™ D = 1.99 ±0.06 . (95) 

Alternatively, it is possible to determine the pion-photon coupling directly. 
Analogous to Eq. (^) , we can write down a gauge-invariant photon-emission 
transition, 

X 7 = ^ie^p,ti{k)q a h P(7 q° . (96) 
m A 

Fixing the coupling g nl in terms of the decay rate 
yields the value 

p^p* = 0.98 ±0.11 . (98) 

Thus, one obtains a factor-of-2 discrepancy between the empirical ampli- 
tude and the VMD prediction, with the VMD value being the larger. Several 
possible explanations for the lack of agreement come to mind. Although 
the radiative branching ratio given in Eq. (^) has reasonably small error 
bars, the signal is based on only one experiment. Alternatively, there may 
be unexpectedly large momentum dependence in the A 2 7rp vertex. Thus, as 
one proceeds from the rho mass-shell (k 2 = m 2 ) to the photon mass-shell 
(k 2 = 0), a 'softening' might occur in the VMD estimate. However, to our 
knowledge there is no previous evidence for such momentum dependence for 
the p extrapolation. 
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